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In the Cold Dark Matter (hereafter CDM) scenario even isolated density peaks contain 
a high fraction of small scale clumps having velocities larger than the average escape 
velocity from the structure. These clumps populate protoclusters, especially in the pe¬ 
ripheral regions, r > Rf (where Rf is the filtering scale). During the cluster collapse and 
the subsequent secondary infall, collapsing or infalling clumps (having v < Uesc) interact 
with the quoted unbound clumps (or high speed clumps, as we also call them) having 
V > Uesc. We study the interaction between these two kind of clumps by means of the 
impulse approximation^ and we find that the collapse of bound clumps is accelerated 
with respect to the homogeneous case (Gunn ^ Gott’s model. Ref. 2). The acceleration 
of the collapse increases with decreasing height of the peak, v. We finally compare the 
acceleration produced by this effect to the slowing down effect produced by the gravi¬ 
tational interaction of the quadrupole moment of the system with the tidal field of the 
matter of the neighboring proto-clusters studied in Del Popolo Sz Gambera^. We find 
that the magnitude of the slowing down effect is larger than the acceleration produced 
by the effect studied in this paper, only in the outskirts of the cluster. We want to stress 
that the one which we study in this paper is also present in an isolated protocluster, 
being produced by the interaction of the collapsing clumps with the unbound substruc¬ 
ture internal to the collapsing clumps itself while that studied in Ref. 3 is produced by 
substructure external to the density peak. 

Keywords: cosmology: theory-large scale structure of Universe - galaxies: formation. 
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1. Introduction 

According to the most promising cosmological scenarios, structure formation is 
traced back to the evolution of primordial density fluctuations. These fluctuations, 
originated from quantum fluctuations^’^’®’^ in an inflationary phase of the early 
Universe, grew up through gravitational instability to a maximum radius r-^. At 
the time corresponding to maximum expansion, perturbations broke away from 
the general expansion and at (5 ~ 1 them began to collapse. Hence the collapse of 
perturbations onto local density maxima of the primordial density field has a key 
role in structure formation and several studies deal with this problem. The problem 
of collapse has been investigated from two points of view, namely: 

i) that of the statistical distribution of the objects formed"^’® 

ii) that of the structure of these objects and its dependence on the statistical prop¬ 

erties of the primordial density field^’®’^®’^^ 12.13,14.15 16,17, is, 19 20,21,22,23 

24 


In the spherical accretion model introduced by Hoyle and Narlikar^® and applied to 
clusters of galaxies by Gunn & Gott^, the matter around the core of the perturbation 
is a homogeneous fluid with zero pressure. If the density inside the perturbation is 
greater than the critical density, it is bound and shall expand to a maximum radius 
I’m: 


I’m = 


5 


( 1 ) 


where ri is the initial radius and 5 is the overdensity inside the radius r. Such matter 
shall collapse in a time 


Tco/2 


TT (I -f l5) 


( 2 ) 


where Hi is the Hubble parameter at the initial time U. 

This model was introduced in order to overcome the problem of the excessively steep 
density profiles, p oc r“^, obtained in numerical experiments of simple gravitational 
collapse. Some authors were able to produce shallower profiles ^p ex r“^, 
through the secondary infall process. Several years later, observational evidences 
for secondary infall in the outskirts of clusters of galaxies has been reported in Ref. 
28, 29. 

Even if this model and in general the original form of SIM (secondary infall model) 
is able to explain better than previous models the structure of clusters of galaxies 
it has some drawbacks, (e.g. it does not predict the right structure of density pro¬ 
files), that can be overcome with some improvements in the original model®®. Two 
noteworthy limits by the Gunn & Gott^ model are the following: 


a) It neglects tidal interaction of the perturbation with the neighboring perturba¬ 
tion; 
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b) It neglects substructure existing in the outskirts of proto-galaxies or proto¬ 
clusters and in the background. 

The effect of the tidal field was studied in a recent paper^. In that paper, we showed 
that the gravitational interaction of the quadrupole momentum of a proto-structure 
with the tidal field of the neighboring proto-structures introduces another potential 
energy term in the equation describing the collapse which acts in the sense to delay 
the collapse. This effect has revealed as one of noteworthy importance in solving 
several of the drawbacks of the CDM model^’^^’^^’^^, namely: 

1) lack of a mechanism originating the bias^; 

2) discrepancies between: 

• al) X-ray temperature function of clusters and observations^^. 

• a2) two-point correlation function of clusters and observations^^; 

• a3) mass function and velocity dispersion function of clusters and 
observations^^; 

• a4) shape of clusters and observations^"*; 

• a5) angular two-point correlation function of galaxies and observations^^. 

For what concerns point ” &)”, we may say that during the last two decades, as¬ 
tronomers have discovered that at least a third of galaxy clusters are not dynami¬ 
cally relaxed systems but contain substructure on scales of the same order as the 
cluster itself^"*This implies that clusters are currently forming or have formed 
recently enough that they have not had time to undergo significant degree of violent 
relaxation and phase mixing. When one studies the structure of the velocity field 
around a peak one finds that at radii r larger than the filtering scale R{, (r > R{), a 
sensitive fraction of the matter of the peak is gravitationally unbound to it^^. This 
is not surprising since a non zero fraction of unbound matter should be expected 
even in the case of a uniform system with a Maxwellian velocity distribution (see 
section 2.2 for a discussion). 

The encounters between collapsing clumps and the unbound high velocity clumps 
influence the collapse of the proto-cluster. We want to study this point and to com¬ 
pare its effect on cluster collapse with that produced by the effects of tidal interac¬ 
tion of the protocluster with the neighboring ones. 

In the following we show how the unbound substructures influences the collapse 
of the perturbation. The plan of this work is the following: in Section 2 we deal 
with the unbound substructure effects. In Section 3 we give our results and finally 
Section 4 is devoted to the conclusions. 

2. Unbound substructure effects 

Our cosmological setting is described by a hierarchical clustering scenario, a CDM 
model (see Ref. 38) in which structure formation goes bottom-up. Local density 
extrema are assumed to be the sites of structure formation^’®. The perturbations 
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grows linearly until the linearly evolved dispersion < 7 ( 2 ;, M) of the density contrast 
is < 1, while it enters in the non-linear phase when a{z,M) = 1. Denoting this 
epoch by Zni one has: 


Cro(M) = 1 -I- Zni 


( 3 ) 


(Liddle & Lyth^®, equation 8.2), where as always the subscript 0 denotes the present 
value of the linearly evolved quantity and cto(M) is related to the power spectrum 
P{k) of density fluctuations by: 

ao(M) = ^ J P{k)k'^W{kRi)dk (4) 


being i?f the filtering scale and W{kRi) the window function: 


W(kRi) = 


3[sin{kR{) — kRfCos{kR{)] 


A subgalactic mass scale collapses at a redshift: 

1 30 

1 + ^nl — ^ 

0 


( 5 ) 

( 6 ) 


(Silk & Stebbins®®), where b is the biasing parameter. In particular a perturbation 
of IO^^Mq collapses at a redshift z ~ 0.02 while perturbations in the range IO^MqA- 
10^Mq collapse almost at the same 2 ; ~ 18 because the mass variance in this region 
varies only of a factor 3 (Ref. 40). A part of this last perturbations have a cross 
section too little for gravitational merging to be destroyed"^®. Therefore a mass 
perturbation > 10 ®Mq (see Ref. 38) can survive until the cluster enters non-linear 
phase at 2 ; ~ 0.02. 

In order to study the dynamical evolution of the proto-cluster, we use the shell 
method®If we consider a shell of matter around the centre of the cluster this 
is made of galaxies, gas and substructure. During the collapse, the shell encounters 
clouds of substructure some of which have high speed. Supposing that the systems 
subject to the encounter (the shell of matter and the high speed clumps) have, 
respectively, median radii ri and r 2 , velocity dispersions of order cti, (J 2 , that at 
the instant of closest approach their centres are separated by P and have relative 
speed V the condition"^^: 


max{ri,r 2 ,b') 
V (J'l - 


i=l,2 


( 7 ) 


ensures that in the course of the encounter the objects constituting the encounter¬ 
ing systems barely move from their initial locations. The result of the encounter is 
to generate perturbations in the velocity and energy of the objects constituting the 
shell (here we are only interested in the effect of the encounter on the shell). In fact 
the effects of high speed encounters on the internal structure of ” stellar” systems 
decrease as the encounter speed increases and consequently high speed encounters 
generate only small perturbations of otherwise steady-state systems. 
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The conditions expressed in equation ( 0 ) is known as impulse approximation (here¬ 
after lA)^ and yields accurate results often even when the condition is not strictly 
satisfied. As stressed by Gnedin & Ostriker"*^, the IA works well in the outer region 
of clusters, wherein galaxies are assumed to move little during the interaction, while 
the lA is not appropriate in the cluster core. 

2.1. Applicability of the I A 

Before going on, we want to analyze the applicability of the lA as expressed by 
equation (Q) to our problem. 

A first rough analysis can be performed remembering that the escape velocity 
in the cluster is fixed by the virial theorem according to equation (see for example 
Ref. 44, (equation (8.3))): 

< >=4.<v^> ( 8 ) 

where Vesc is the escape velocity and < > is the 3-D dispersion velocity, and 

< >= 3al (9) 

being cti the velocity dispersion of the objects constituting the shell, and then 

■Uesc = '/T2(7i (10) 

The relative velocity, V, is the sum of the infall velocity of the shell and the velocity 
of clumps moving outwards: 

R — r^clumps 4“ "^infallingshell (H) 

Imposing that the clumps move at speed larger than Vesc, we get: 

^clumps ^ ^esc ~ '\/l2o'l = 2V^Cri (f^) 

The infalling shells are marginally bound, hence they travel at a velocity close to 
the local escape velocity 

^infalling shell ^ ^esc — VUai (13) 

This last statement can be easily proven using for example an Hernquist model for 
the cluster. In general, conservation of energy yields: 

^infallingshell ^esc 4” (1'^) 

where $ is the (negative) potential and is the maximum (apocentric) radius of 
the infalling particle. 

For example, for a Hernquist cluster, one has 

$(r) = —GM/(r 4-a) (15) 

and 


^^hifaiiing shell (?-) = 2[$(ra)-4>(r)] = (r) [1-(r-ko)/(ra-ko)] = v^^^(r)(ra-r)/(ra+a)(16) 
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For a Hernquist model, a is a few tenths of the virial radius. The infalling 
particles are between the virial and turnaround radii, i.e. between rvir and roughly 
2 Tvir- So Tq >> o and 

^infallingshell(^) ^esc(^)(^ ^/^a) 

which means that r’infaiiing shell is a little smaller than v^sc- 

Using the previous results (equation (|l^and ), we have: 

U — riclumps “f 'llinfalling shell — 2\/3(7i + "s/120'i (1^) 

If z = 1, in Eq. (|^), we are dealing with the shell then: 

1) U = 7ai > CTi 

If z = 2, we are dealing with the clump: a 2 — Ikm/s << CTi (see Ref. 39) and then: 

2) U = 7cti >> {72. 

Both the relations 1) and 2) satisfy the condition: 

V > 7ai (19) 

value that according to Ref. 44 (see page 440) is considered good to use the lA. 
We want to remark that the previous analysis does not take account of the reduction 
of escape velocity with radius: given the potential of the sistem, 4>(r), the escape 
velocity is Vesc = \/2|$(r)|. This means that going outwards the lA is more strictly 
satisfied. 


2.2. Energy change in the shell 

The energy change in the shell, due to an encounter with high velocity clumps, can 
be computed as follows. 

Suppose to have a perturber (a high speed clump of substructure, member of the 
cluster) characterized by a speed larger than the escape speed from the cluster. 
This high velocity clump, that we stress is inside and part of the cluster, moves 
away from it. We use cylindrical coordinates (R,z) such that the z-axis coincides 
with the perturber’s trajectory given by: 

i? = 0, z = A{t) = Vt (20) 


and that i? = 0, z = 0 is the centre of the perturbed system (shell of matter). Note 
that we are interested only in clumps not contained in the core and moving from 
the centre in the outward direction, so that they do not go through the cluster 
centre. The choice of not taking into account the effect of clumps in the core area is 
due to Gnedin & Ostriker"^^ observation that the lA works well in the outer region 
of clusters, wherein galaxies are assumed to move little during the interaction. The 
velocity increment Av of a test particle at the position (R,z) is given by: 


Av^ = 


d<^R R 


dd) dA 
dr r 


0 


0 


( 21 ) 
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where 

r = [(A-z)2 + i?2]i/2 (22) 

and <i> is the potential of the perturber. Assuming that the perturber is an isothermal 
sphere it can be easily shown that: 

Awr = —TTa^/V (23) 


where CT 2 is the perturber internal velocity dispersion. 

The energy, gained by the shell of matter during the encounter with a high velocity 
clump, is given by: 

poo 

AE= [Av{R)fE{R)RdR (24) 

Jo 


(see Ref. 44) where S(i?) is the surface density of the perturbed system. 

We assume that the perturbed system (shell) is characterized by an average 
density given by: 


p(n,t) 


3M 

47ra^(n, t)rf 


(25) 


where, according to Ref. 2 notation, a{ri,t) is the expansion factor of the shell at 
time t, Ti its initial radius and M its mass. We are studying the collapse of a shell of 
matter then and ri coincide with the turn-around epoch and radius of the shell, 
respectively. Assuming that no shell crossing occurs, the total mass inside each shell 
stays constant and then: 


p(n,t) = 


a3(ri,t) 


The surface density is given by: 


E = 


Si 


a(n,t)2 


(26) 


(27) 


At this point we are ready to calculate the energy change in the shell by using 
equations dl), (H), (izl) and remembering that r = ria{ri,t): 


AE = 7r2-|^Eirf(logai„ax - logOmin) 


(28) 


where Umax is the maximum value of the expansion parameter, corresponding to 
the turn-around epoch, and Omin the minimum value of the expansion parameter 
corresponding to the time of full collapse of the perturbation. During infall a shell of 
matter experiences several encounters with all the unbound clumps located between 
its initial position and the final one. Each one of the encounters produces an increase 
of the energy of the shell. In order to calculate the total energy, A At, acquired 
during infall we must calculate the number of unbound clumps ’’internal” to its 
orbit. The number of perturbers can be calculated observing that: 
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aa) they are peaks of the density field having dimension much less than that of 
the cluster but large enough to survive till the cluster enters the non-linear 
phase; 

bb) they must satisfy the lA condition. 


Clumps in the SCDM scenario are not distributed randomly but are clustered and 
so tend to be more bound to their nearby neighbors. This may make them more 
subject to mergers in some cases soon after they form, but may also prevent mergers 
by placing them in a high velocity dispersion environment. For survival, we must 
require a large density contrast to have developed. If the smaller scale fluctuations 
have a density contrast of ~ 100, they should be able to survive tidal disruption 
or dynamical friction drag when its environment undergoes non-linear collapse. 
This translates into a minimum mass scale of ~ IO^Mq (see Ref. 39). In order 
to determine the number density and the characteristics of these peaks we follow 
the same choice made by Antonuccio-Delogu & Colafrancesco'^'* and Del Popolo & 
Gambera^^. We restrict ourselves to consider a subset of peaks of the density field, 
having the central height v larger than a threshold v^- The threshold is chosen by 
satisfying the condition its radius rpk is smaller than the average peak separation, 
dav = > Vc): 

rpk(i^ > Vc) < 0.1 • > Vc) (29) 


where the radius of a peak rpk is given by: 

rpk = ^ 

(see Ref. 45), where: 


1 + z/f7o 7^ + (0.9z/)^-^J 

7 = all ao(T2 


(30) 


(31) 


and here Ui is the i-th momentum of the variance of the fluctuation field: 

erf = ^ J P(fc)fc20+l)rffc (32) 

where P{k) is the filtered power spectrum. The power spectrum that we adopt is 
P{k) = AkT'^{k) with the transfer function T{k) given in Ref. 7 (equation (G3)): 

T{k) = + (16.1g)2 -h (5.469)3 -h {6.71)^^^^ (33) 

where A is the normalizing constant and q = Here 0 = per/(l-68p.y) 

represents the ratio of the energy density in relativistic particles to that in photons 
(0 = 1 corresponds to photons and three flavors of relativistic neutrinos). For 
the previous GDM transfer function, normalized as as = 0.63 and filtering radius 
R{ = 30Kpc, (corresponding to a Gaussian filtered mass M ~ 10®Mq) one has 
7 = 0.48, i?* = = HKpc and then Vc = 1.6. Remembering that the number 



F(v>V 
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Fig. 1. Spatial variation of F{v > t;esc) for the CDM spectrum, introduced in Section 2, smoothed 
on a scale = bh~^ Mpc and for u = 1.2. 


density of peaks, ria, having v > is given by: 
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Fig. 2. Energy acquired during encounters of shells of matter infalling from a distance r from the 
centre of the proto-cluster. The dashed line represents AEx for the CDM spectrum, introduced 
in Section 2, smoothed on a scale Rf = Mpc and for = 1.2. The solid line is the same but 

now V = 2 


where 

^pk = 7:7^^G'(7,7J^)exp(-i/V2) 


(35) 
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being G( 7 , 7 ^) defined in Ref. 7 (equation (4.4)) together with and 7 , we obtain 
ria — 135Mpc“^ The average mass of the peaks with v > i/c is given by: 


where 


i(l/ > t'c) 


pOC) 

/ TOpeak^pk(l")di^ - 7 X 10 ®Mq 

J I/c 


_ _ 47r ^3 2 ®/^ 

mpeak - 


(36) 


(37) 


(see Ref. 46). This mass satisfies both condition ” aa) ” and Silk & Stebbins®® 
prescription for surviving clumps. 

The fraction of this subpeaks satisfying the lA, and consequently satisfying condi¬ 
tion ” bb)”, can be calculated using the conditional probability distribution /pk(v|i/) 
giving the peculiar velocity around a peak of a gaussian density field having central 
height V. 


/pk(v|:^) = 




• e 


-EL 


3 




(see Ref. 37) (<> denotes ensemble averages), where 

v(vx) 


V2 . 

a =- erjc 

TT 


{yx) = 7 , (Ref. 7) 


(vv,) = -- 


HM-^n 


^(r) is the autocorrelation function and 

a; =-- 

O’ 2|0 


L -y/Ar , 

dxx^^{x) 


(38) 


(39) 


(40) 


(41) 


is the central curvature of the peak. Ai, A 2 , A 3 are three of the four eigenvalues 
of the covariance matrix M described in Antonuccio-Delogu & Colafrancesco®^. 
These three eigenvalues represent the components of the peculiar velocity disper¬ 
sions along directions parallel (Ai) and perpendicular (A 2 ,A 3 ) with respect to the 
radial direction of a density peak. The fraction of peaks having v > Vesc can be 
obtained by means of the conditional probability distribution /pk(w|j^): 


F{v > Vesc) 


fo°° fpk(vli^)d^v 


(42) 


The fraction of peaks satisfying the condition F(v > Vesc) can be easily calcu¬ 
lated to be F{v > Vesc) = 0.00738 for a Maxwellian velocity distribution, f{v) oc 
exp(—^^), characterized by a uniform velocity dispersion a^. In the case of a top- 
hat perturbation of radius Rf = 1.5/i-^Mpc, F{v > Vesc) can be as large as 0.3 for 
a typical velocity dispersion, and even larger for larger values of R{. 
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Using the conditional probability distribution /pk(v|i/), in the case of a = 1.2 
peak we get F(v > Vesc) — 0.09 — 0.13 within iif = 5/i“^Mpc. 

Q In Fig. 1, we plot the spatial variation of F(v > Vesc) for the CDM spectrum, 
introduced in Section 2, smoothed on a scale Ri = 5h~^ Mpc and for i/ = 1.2. 
In any case in the following, we’ll make a conservative calculation using the lower 
estimates for F{v > Vesc) — 0.09. 

Finally, the total number of subpeaks verifying the above conditions ” aa) ”, ” 
bb) ” is given by: 

47r o 

Ntot = • na • F{v > Vesc) (43) 

while the total energy, AE^, injected in the shell by the impulsive encounters with 
all high speed subpeaks is given by: 


AE'y — Ntot • 


(44) 


The energy for unit mass of a collapsing shell is a modified version of that given by 
Ref. 47 (equation( 19.9)) and can be written similarly to Ref. 48 in the form: 




GM 

r 


AE^ 


(45) 


where M is the mass of the shell of matter. 

Integrating equation (^), using equations (|^, (^) and remembering that pi = 
Pci(l + iJ) where pci = with pd and H{ respectively the critical mass density 
and the Hubble constant at the time ti, the time of collapse is given by: 


A = 2 


da 


H? 


l+(5 _ l+(5 


4G 


Hf{l+S)rf 


AUt 


Integrating analytically equation ([461), we find the solution: 

TT 1 

ic/2~ 




(s + ^4^5 aeAj 


•ill 


(46) 


(47) 


having assumed <5 << 1. In the case AAt = 0 the solution of this equation is that 
given In Ref. 2: 


ico/2 


TT 1 

-H.jn 


(48) 


^It is important to note that the use of the peak formalism is not fundamental for what concerns 
the results of this paper. Similar results are obtained even using a top-hat perturbation for which, 
has previously remarked, F{v > Uesc) can be as large as 0.3. 
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Fig. 3. The time of collapse of a shell of matter around a peak having i' = 1.2 in a CDM model 
taking account of the effect of the unbound high velocity clumps (solid line) and for the Gunn Sz 
Gott^ model (dashed line), to is Hubble time. 


3. Results and discussion 

The results found by means of the model introduced in the previous section are 
shown in the Fig. 1-5. 




February 1, 2008 21:8 WSPC/Guidelines-IJMPA ijmpa 


14 A. Del Popolo, M. Gamnera, E.N. Ercan 


As previously reported, Fig. 1 shows the spatial variation of F(v > Uesc) for 
the CDM spectrum, introduced in Section 2, smoothed on a scale Rf = Mpc 
and for v = 1.2. The fraction of clumps with v > r’esc within Rf = 5/i“^Mpc is 
F{v > Vesc) — 0.09 — 0.13, and increases with increasing radius. This behavior is 
due to the fact that the potential of the structure becomes shallower going towards 
its outskirts, and consequently Vesc is substantially lower. In Fig. 2 we show the 
total energy injected in the shell, because of the impulsive encounters of collapsing 
shells with the high speed clumps, in function of the initial radius of the shell. The 
dashed line represents AEt for a CDM model with spectrum given by equation ( ^3| ) 
for V = 1.2. As shown, shells collapsing from larger distances from the centre of 
the cluster acquire more energy. This is due to two different reasons: firstly the 
energy, AA, acquired in an encounter with a single clump increases with increasing 
turn-around radius of the shell of matter. Secondly, AAt, is proportional to the 


number of unbound clumps that increases with radius (see equation (43)) and then 
a shell falling from a larger distance encounters more unbound clumps. The solid 
line represents A Ax for z/ = 2. The value of AAx is smaller with respect to the 
previous case. The reason is due to the fact that the fraction of the unbound peaks 
decreases for higher v peaks because the potential of these protostructures is deeper 
and hence the value of Vesc is systematically higher. 

In Fig. 3 we show the collapse time versus the distance from the centre of the cluster. 
The figure plots the collapse time from 1 /i“^Mpc on, because of the limits of the 
lA, stressed by Ref. 42, in the central part of the cluster. The solid line represents, 
tc, for = 1.2 and taking account of the unbound clumps effect, while the dashed 
line is the collapse time, tco, given by the model in Ref. 2. This figure shows that 
the collapse time is accelerated, with respect to the model in ref. 2, because of the 
energy acquired by the shells during encounters with unbound clumps. In Fig. 4 
we show tc, for v = 2. As in the previous figure the solid line represents R taking 
account of the unbound clumps effect, while the dashed line represents tcO given 
by the model in Ref. 2. Comparing Fig. 3 and Fig. 4 we see that the effect of the 
unbound peaks on the acceleration of the collapse time is larger in the case v = 1.2. 
This is in agreement to what observed for AE^. 

Finally, in Fig. 5 we make a comparison between the acceleration in collapse time 
due to unbound clumps to the slowing down effect produced by the gravitational 
interaction of the quadrupole moment of the system with the tidal field of the 
matter of the neighboring proto-clusters, described in Ref. 3. 

As shown in Fig. 5, the magnitude of the two effects is comparable at lh“^Mpc, 
but the accelerating effect produced by the high speed clumps decreases with in¬ 
creasing distance from the cluster centre, with respect to the slowing down effect 
produced by the tidal field. 

The mechanism considered in this paper is different from that described in Ref. 3 
for two reasons: 


• the delay of the collapse produced by unbound clumps is due to encounters 
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Fig. 4. Same as Fig. 3 but now u = 2 


of infalling shells of matter and unbound clumps. The effect described in 
Ref. 3 is produced by tidal interaction of the matter of a given protocluster 
with that of the neighboring ones; 

• the substructure responsible for collapse delay in this paper is that con¬ 
tained inside the infalling shell while in Ref. 3 it is external to the shell. 
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Log r (h ^ Mpc) 


Fig. 5. Comparison of the collapse time Tc^unbound)? obtained taking account unbound high 
velocity clumps, with Tc(tidal)? obtained taking account of the tidal interaction of the protocluster 
with the neighboring ones. In this case v = 2 


The two effects have a different physical origin, and act in the opposite direction: 
unbound clumps accelerate the collapse, the effect of the tidal field is to delay 
the shell collapse. It is obvious that during the infall of a shell both effects are 
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present with the final result that at small radii (~ l/i“^Mpc) they counterbalance 
while at larger radii tidal field is dominant producing a slowing down of the shell 
collapse. We want to stress that the idea that high velocity clumps can influence 
the collapse and secondary infall of matter in protoclusters was introduced (but not 
developed) in a paper by Antonuccio & Colafrancesco^^. In the conclusion of that 
paper, regarding the velocity correlation tensor and unbound substructure around 
density peaks, they noted that high velocity clumps can exchange momentum with 
the collapsing or infalling clumps and generate a drag force onto the collapsing 
material with a consequent delay of the collapse. According to their paper the drag 
force is due to dynamical friction. Unfortunately this last considerations are not 
quite right because: 

1) the effects of the high velocity clumps on the collapse shell of matter cannot be 
calculated using the dynamical friction approximation since a necessary condition 
to apply this last is that the encounter of two systems is characterized by U ~ ftip 
(where Utip is the internal velocity of the system) (see Ref. 44). To be more precise, 
according to Ref. 49, a body of mass M moving through a homogeneous and infinite 
distribution of bodies of mass m, that we call field particles, is subject to the 
deceleration: 


dVM 

dt 


— \nAG^m{M + m) 


/J'“ f{vm)vl,dvn 


-Vm 




(49) 


where In A is the Coulomb logarithm and /(fm) is the phase-space number density 
field particles. As clearly stated by equation (^ ) only field particles moving slower 
than M contribute to the force. In our case, the role of the field particles is played 
by the high speed clumps that have speeds much larger than those characterizing 
the motion of the shell particles (when lA is satisfied) or at least a bit larger. 

2) The effects of the high velocity clumps is to accelerate the shell collapse and not 
to delay it. 

We note that the erroneous conclusion of Antonuccio & Colafrancesco^^ paper is not 
connected with their analysis of the conditional probability distribution fp\^{'v\v) 
developed in the paper, that we used to calculate the number of subpeaks satisfying 
the lA. The erroneous conclusion is due to the way they thought to calculate the 
interaction between the high velocity clumps and the infalling shells. 


4. Conclusions 

In this paper we have studied the role of unbound high velocity clumps on the 
collapse of density peaks in a SCDM model (Uq = 1, h = 0.5, n = 1). We have 
shown that: 

1) the encounters of collapsing shells with unbound clumps produce an acceleration 
in the collapse of density peaks. The effect is larger for peaks having a lower value 
of v (see Fig. 2-^3). 

2) A comparison between the collapse acceleration due to encounters of shells with 
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unbound clumps and the slowing down effect produced by tidal interaction between 
a protocluster and the neighboring ones (see Ref. 3) has shown that these two effects 
act in opposite direction, namely the first accelerate the shell collapse, while the 
second delay the protocluster collapse, especially in its outskirts. The magnitude of 
the two effects is comparable at l/i“^Mpc but the effect of the tidal field dominates 
at larger radii. 

3) The acceleration in the collapse is due to the additive action of several high speed 
clumps with the matter of the infalling shells. 
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